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Abstract. For the first time a PbW04 crystal was grown using ancient Roman lead and it was run as a 
cryogenic detector. Thanks to the simultaneous and independent read-out of heat and scintillation light, 
the detector was able to discriminate /J/7 interactions with respect to a particles down to low energies. 
New more stringent limits on the a decays of the lead isotopes are presented. In particular a limit of 
T 1/2 > 1.4-10 20 y at a 90% C.L. was evaluated for the a decay of 204 Pb to 200 Hg. 

PACS. 23.60.+e a decay - 29.40.Mc Scintillation detectors - 07.20.Mc Cryogenics; low-temperature de- 
tectors 



1 Introduction 

After the first observation of 209 Bi [1J a decay, and the 
recent measurements of the half-life value including the 
transition to ground and to the first excited state [2] , lead 
is considered to be the heaviest stable element. However 
a decay in lead is energetically allowed for all the four 
naturally occurring isotopes. The Q-values of the decays 
(Q Q ) and natural abundances (S) are listed in Tab. m 
as well as calculations of half-lives based on the cluster 
models (3j|4| , phenomenological fission theory of alpha de- 
cay [5] , semiclassical WKB approximation [6 , and micro- 
scopic approach [7]. Given the values in Tab. [T] there is 
no feasible perspective to observe the a decay of 206 Pb, 
207 Pb and 208 Pb, but the measurements described in the 
present work might test the reliability of the various nu- 
clear models. The same perspective holds also for 204 Pb , 
but in addition, this isotope plays also an important role 
in lead geochronology, as a standard reference [8]. 

In 1958 there was a claim [9, of a possible a decay 
of 204 Pb. The nuclear emulsion technique with enriched 
(27.0%) 204 Pb was used for this experiment. A peak was 
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found between 8 ^m and 9 fj,m in the emulsions |10| . cor- 
responding to an a particle of « 2.6 MeV. This result was 
excluded by the mass data on 204 Pb and 200 Hg and 
the same measurement was then translated to the present 
quoted limit: T 1/2 > 1.4-10 17 y. 

The need of increasing the experimental sensitivity to 
rare a decays lead to the employment of different and new 
experimental techniques such as ionizing chambers |12j . 
crystal scintillators [13] and liquid scintillators |14| . 

Recently, measurements of a decays of 209 Bi [TJ|5] and 
180 W P~5|, with half-lives of 1.9-10 19 y and 1.8-10 18 y re- 
spectively, demonstrated the high sensitivity for the dis- 
covery of rare nuclear processes that can be achieved with 
scintillating bolometers. The main advantage of this tech- 
nology is the wide choice of detector materials that allows 
to investigate isotopes that are not easily measurable with 
conventional detectors. Moreover, the simultaneous read- 
out of light and heat signals results in a powerful tool for 
background identification, thanks to the different amount 
of light emitted by a and /J/7 particles of the same en- 
ergy. This is fundamental in case the expected energy of 
an a decay lies in the environmental /J/7 energy region 
(dominated by the 2615 keV 7 line of 208 T1). 
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Table 1. a transitions in Pb isotopes to the ground state of 
daughter nuclei. Q a values are taken from [16] . natural isotopic 
abundances (<5) from [17] . The theoretical expected half-lives 
are evaluated using theoretical models given in |3]|4][5]|6][7] . 



Isotope 


8 


Q Q 


rptheor 




[%] 


[keV] 


[y] 


204p b 


1.4 (1) 


1969.5 (12) 


2.3- 10 35 -^ 1.2- 10 37 


206p b 


24.1 (1) 


1135.5 (11) 


1.8- 10 65 -^ 6.7- 10 68 


207p b 


22.1 (1) 


392.3 (13) 


3.6-10 152 ™ 3.4- 10 189 


208p b 


52.4 (1) 


516.9 (13) 


1.2-10 124 -f 7.4- 10 132 



2 PbW0 4 crystal from ancient Roman lead 

When coming to low background measurements, the main 
problem of lead is the high activity in 210 Pb. The radioac- 
tivity of "commercial" lead can be of the order of few tens 
of Bq/kg, for ore-selected samples [T5] > while on other sam- 
ples it can reach an activity of thousands of Bq/kg [TU] . 

Besides the excellent performances that characterize 
thermal detectors, a limiting factor arises from the rela- 
tively slow time response of these devices. Depending on 
the technology of the temperature sensors, as well as on 
the dimension of the absorber, the time development of 
the thermal pulse can range between few /jjs up to few 
seconds. This last point has to be taken into account in 
the case of crystal compounds based on lead. 
210 Pb is a decay product of the 238 U decay chain, which 
is present in all rocks and ores. During the melting of the 
ore, 210 Pb concentrates in the lead-metal, while all the 
other long-lived radioactive nuclides of the 238 U chain ( 
234 U, 230 Th and 226 Ra) are extracted from the slag since 
they are chemically different from Pb. Since the 210 Pb 
half-life is 22.3 y, its activity should be extremely small in 
ancient samples. For this reason ancient Roman lead [20] 
was used to grow our crystal. This lead, in form of ingots, 
was recovered from the wreck of a Roman ship sunk near 
the coast of Sardinia in the Mediterranean Sea. The 210 Pb 
content of this Roman lead was measured to be less than 
4 mBq/kg. We decided to grow a standard lead containing 
scintillator like a PbWC>4 crystal [21]. A small sample of 
this lead was sent to the Preciosa company in Hungary, to 
grow an undoped PbWC>4 crystal. The request to have a 
clean and undoped crystal arises from two different con- 
siderations: 

— in a bolometer, the thermal signal strongly depends 
upon the crystal purity; 

— the scintillation yield of doped crystals normally de- 
creases at very low temperatures, while for pure crys- 
tals it increases |22| • 

The PbWC>4 crystal used in this work has dimensions 
of 3.0x3.0x6.1 cm 3 and a total mass of 454.1 g. A tiny 
splint (~ 50 mg) was removed from the crystal and ana- 
lyzed through ICP Mass Spectroscopy in order to evaluate 
the isotopic abundances of the four lead isotopes. The re- 
sults are shown in Tab. [2] 



Table 2. Lead isotopic abundaces for the PbWCU crystal, 
evaluated with ICP-MS measurements. The values are in good 
agreement with the natural abundances quoted in Tab. [I] 



204p b 206 Pb 207 Pb 208 Pb 

1.34±0.02 % 25.10±0.25 % 21.1±0.30 % 52.4±0.50 % 




Fig. 1. Set-up of the detector. Two Cu columns, one fixing 
the S-shaped PTFE and one fixing the LD Cu frame are not 
visible due to the chosen cross section. 

3 Experimental Set-up 

The scheme of the set-up is shown in Fig. [I] The PbW0 4 
crystal is held by means of four S-shaped PTFE supports 
fixed to four cylindrical Cu columns. It is surrounded 
(with no direct contact) by a plastic reflecting foil (3M 
VM2002). The Light Detector (LD) [27| is constituted 
by a 36 mm diameter, 1 mm thick pure Ge crystal ab- 
sorber working as a bolometer: it is heated up by the 
absorbed photons and the temperature variation is pro- 
portional to the scintillation signal. The Ge wafer is held 
by two custom-shaped PTFE holders embedded in a cylin- 
dric Cu frame. The frame is held above the PbW04 crystal 
using two Cu columns. 

The temperature sensor of the PbW04 crystal is a 
3x3x1 mm 3 Neutron Transmutation Doped (NTD) ger- 
manium thermistor, the same used in the Cuoricino ex- 
periment [28J. It is thermally coupled to the crystal via 
9 glue spots of « 0.6 mm diameter and w 50 /im height. 
The temperature sensor of the LD has a smaller volume 
(one third of the PbW04's one) in order to decrease its 
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heat capacity, increasing therefore its thermal signal. The 
PbWC>4 detector was run in an Oxford 200 3 He / 4 He dilu- 
tion cryostat deep underground in the Laboratori Nazion- 
ali del Gran Sasso. The details of the electronics and of 
the cryogenic facility can be found elsewhere [2"9"II5U] , 

The heat and light pulses, produced by interacting par- 
ticles in the absorber (transduced in a voltage pulse by the 
NTD thermistors) are amplified and fed into an 18-bit NI- 
6284 PXI ADC unit. 

The trigger is software generated on each bolometer 
and when it fires waveforms 1 s long, sampled at 1 kHz, are 
saved on disk. Moreover, when the trigger of the PbWC>4 
crystal fires, the corresponding waveform from the LD is 
recorded, irrespective of its trigger. 

The amplitude and the shape of the voltage pulse is 
then determined by the off-line analysis. 

The heat axis is energy-calibrated attributing to each 
identified peak (due to an external 7 source and to internal 
a contaminations) the nominal energy of the line. Two 
independent functions were used for calibrating the a and 
(8/7 bands. 

The dependency of amplitude from energy is parame- 
terized with a second order polynomial of the heat pulse 
amplitude. 

The energy calibration of the LD is obtained by means 
of a weak 55 Fe source placed close to the Ge wafer that il- 
luminates homogeneously the face opposed to the PbW04 
crystal. 



4 Data Analysis 



Energy [keV ] 

1000 2000 3000 4000 5000 6000 7000 



> 




Energy [keV ] 



Fig. 2. Light vs Heat scatter plot corresponding to 586 h of 
measurement. The horizontal axis reported on the bottom of 
the plot corresponds to an a energy calibration, the one on the 
top a /3/y calibration. 

interest (see Tab. [IJ in order to allow a thorough energy 
calibration. 

One of the constituent of the crystal is tungsten, the 
natural a decay of 180 W at 2516 keV [15] is taken into 
account for the energy calibration. On the other side, the 
LD is calibrated using the 55 Fe X-ray doublet: 5.9 keV and 
6.2 keV. 



Our detector was operated for a total live time of 586 hours 
for a background measurement. In Fig. [2] we present the 
obtained Light vs Heat scatter plot, with the two different 
energy calibrations. 

To maximize the signal-to-noise ratio, the pulse ampli- 
tude is estimated by means of the Optimum Filter tech- 
nique (OF) [31,32 . The filter transfer function is built 
from the ideal signal shape s(t) and the noise power spec- 
trum N(ui). The s(t) is estimated by averaging a large 
number of triggered pulses (so that stochastic noise su- 
perimposed to each pulse averages to zero) while N(lj) 
is computed averaging the power spectra of randomly ac- 
quired waveforms where no pulse was found. The ampli- 
tude of a signal is estimated as the maximum of the filtered 
pulse. The amplitude of the scintillation light signals, in- 
stead, is evaluated from the filtered waveforms at a fixed 
time delay with respect to the PbW04 bolometer, as de- 
scribed in detail in 133 . 



The final spectrum is composed of events which sur- 
vived two different types of data selection global and event- 
based requirements. Global requirements arc applied fol- 
lowing criteria decided a priori on the detector perfor- 
mances (e.g. excessive noise level). They identify bad time 
intervals that need to be discarded. Event-based require- 
ments include: pile-up rejection and pulse-shape selection. 
The presence of a pile-up prevents the OF algorithm from 
providing a correct evaluation of the pulse amplitude. The 
pulse-shape analysis is used to reject non-physical events 
(e.g. electronic spikes). The pulse-shape parameters used 
for selecting the events are the rise time and decay time 
of the OF-filtered waveform and the mean quadratic de- 
viation of raw signals from the average detector response. 



4.1 Efficiency of event-based cuts 



After the application of the OF, signal amplitudes are 
corrected for temperature and gain instabilities of the set- 
up. The PbW04 bolometer is calibrated with the most 
intense a peaks from the internal 238 U and 232 Th contam- 
inations, and with an external 232 Th and 40 K 7 source, see 
Fig. [2] Since we are investigating low energy a decays, the 
internal a contaminations do extend over a large energy 
range, but they are not low enough to reach our region of 



Even if the crystal has been grown using ancient lead, 
the main background source is still 210 Pb due to an in- 
ternal cr ysta l recontamination, as it will be explained in 
the Sect. 4.4 Because of the intense 210 Po contamination, 



there is a significant loss of efficiency due to pile-up event 
rejection (e p u e -up)- The efficiency is estimated as in |34| : 



1 



P, 



pile— up 



-r-T 



(i) 
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Table 3. Efficiencies for a events that passed events-based 
cuts: pile-up, pulse-shape and total. 



Nuclide 


Q- value 


^pile — up 


(-PS 


tTOT 




25f6 keV 


0.94±0.01 


0.88±0.02 


0.83±0.02 




4082 keV 


0.94±0.01 


0.87±0.01 


0.81±0.02 


230 rp^ 


4770 keV 


0.94±0.01 


0.85±0.01 


0.80±0.02 


222 Rn 


5590 keV 


0.94±0.01 


0.87±0.01 


0.80±0.02 


218p 


61 f 5 keV 


0.94±0.01 


0.86±0.01 


0.80±0.02 


210p o 


5407 keV 


0.94±0.01 


0.85±0.01 


0.79±0.02 



Table 4. Technical details for the PbWCU crystal and for the 
LD. Signal represents the absolute voltage drop across the ther- 
mistor for a unitary energy deposition, tr and td are the rise 
and decay time, respectively. 



Crystal 


Signal 


FWHM(, ose 




TD 




[/xV/MeV] 


[keV] 


[ms] 


[ms] 


PbW0 4 


153 


1.51 


13 


35 


LD 


1553 


0.28 


3 


20 



where P p u e -up is the pile-up probability, r is the rate of 
the 210 Po events that passed the global cuts (r = 66 mHz), 
T the time window containing an event during which the 
occurrence of another one would be considered pile-up. 
The acquired waveform are 1 s long, so T is set at 1 s. 
The pulse shape cut efficiency (eps) is estimated at dif- 
ferent energies for various a lines, by a simultaneous fit 
on both spectra of accepted and rejected events |28| . In 
Tab. 3] we report the values for the efficiencies of event- 
based cuts. 



4.3 Light Yields for a and [3/j particles 

As described in Sect. [T] particle discrimination is the key 
point for all the experiments searching for rare decays in 
overwhelming background. 

In Fig. [2] we can distinguish the fi/j and a regions 
which give rise to very clear separate distributions. In the 
upper band, ascribed to /3/j events, the 2615 keV 208 T1 
7-line is clearly visible. The lower band, populated by a 
decays, shows the continuum of events, all the way down 
to few keV, induced by the degraded a of 210 Po, as well 
as 238 U and 232 Th contaminations. 

The Light Yield (LY), defined as the ratio between 
the measured light (in keV) and the nominal energy of an 
event (in MeV), was measured for the various 7- lines and 
it was found to have a constant value: 



LY &h = 1.78 ± 0.01 keV/MeV 



(2) 



On the contrary, the LY a for a decays has a value which is 
correlated to the nominal energy of the interacting parti- 
cle. Fitting the a band with a second order polynomial, we 
can estimate the LY a for the various lines and we obtain 
the following relation: 

LY a (E a ) = (0.28±0.01) + (2.93±0.14)-10" 5 -£; Q keV/MeV 

(3) 

The observed behaviour reflects the fact that an a parti- 
cle produces less scintillation light compared to an electron 
because of its larger dE/dx, which can induce saturation 
effects in the scintillator (the "Birks law" [22])- For the a 
band, the energy dependence of the LY is a consequence of 
the energy dependence of the stopping power. Electrons, 
due to their low stopping power, do not suffer of saturation 
effects and their LY is, consequently, energy independent. 



4.2 Detector performances 



The detector performances are reported in Tab. [4] The 
baseline resolution, FWHMb ase , is governed by the noise 
fluctuations at the filtered output, and does not depend 
on the absolute pulse amplitude. The rise (tr) and decay 
times (td) of the pulses are computed as the time differ- 
ence between the 10% and the 90% of the leading edge, 
and the time difference between the 90% and 30% of the 
trailing edge, respectively. 

The FWHM energy resolution of the crystal ranges 
from 15.2±0.5 keV at 583 keV to 15.7±0.3 keV at 2615 keV 
for the 7 de-excitation of 208 T1. While for the a region it 
goes from 13.8±5.4 keV at 2516 keV ( 180 W) to 52.9±1.9 keV 
at 6115 keV ( 218 Po), both lines refer to a + recoil energies 
of the nucleus. 

The computed FHWM for the 55 Fe on the light detec- 
tor is 379±8 eV. 



The scintillation Quenching Factor (QF) is defined as 
the ratio of the scintillation yield of an interacting parti- 
cle (a, neutron, nucleus) with respect to the LY of a /3/j 
event at the same energy. Since the LY of a events varies 
with the energy of the interacting particle then the QF 
for these type of events does not have a constant value. 
Using the eqs. [2]and|3]we obtain a value of 0.20±0.01 for 
the QF a ( 180 W). This value is in agreement with [25] . 

In our measurement, we observe that two different 
functions are needed for the calibration of the two bands ??, 
as it is shown in Fig. [2] In fact, if we apply the a calibra- 
tion to the (8/7 band, we obtain a large mis-calibration, 
about 80% for the 208 T1 7 line at 2615 keV. The expla- 
nation for this behavior relies in the mechanism of energy 
conversion |25| . In fact, the energy deposited in the abosr- 
ber can be converted into: 

— scintillation light, due to the excitation of luminescence 
centres; 

— thermal phonons; 

— "blind" channels, not detectable. 



J.W. Beeman et al: New experimental limits on the a decays of lead isotopes 



5 




Table 5. Evaluated internal radioactive contaminations for 
the PbW0 4 crystal. Limits are at 90% C.L. 



2500 3000 3500 4000 4500 5000 5500 6000 

Energy [keV ] 

Fig. 3. Energy spectrum of the a region, acquired over 586 h 
of measurement. 



Unfortunately, we are not able to measure the energy par- 
tition in the different channels, since we miss information 
on the light collection efficiency of the LD, and obviously 
we are not sensitive to the "blind" channel. 



4.4 Internal contaminations 



Even if the starting materials for the crystal growth had a 
very high radiopurity level (Sect.[2| during the crystalliza- 
tion process is possible to unintentionally introduce some 
radioactive contaminants. These have been evaluated by 
analyzing the a region of the spectrum, due to the more 
favorable signal-to-noise ratio. 

As already pointed out, the background in this region 
is essentially due to a particles, which produce peaks and 
a continuum of events generated by degraded a's from the 
surfaces of the crystal and the surrounding materials [36 . 
This is mainly due to an intense 210 Po contamination of 
the crystal as shown in Fig. [3] 

By means of Monte Carlo simulations, we estimate the 
background induced by a uniform 210 Po contamination 
throughout the crystal. We obtain a flat background from 
zero energy up to the 210 Po Q-value, the counting rate in 
this region is 0.02 counts /keV. 

We observe that the PbWC>4 crystal is contaminated 
in 230 Th and 210 Po, from the 238 U decay chain, 232 Th 
(the whole 232 Th decay chain) and 147 Sm, see Tab. [5] For 
the 232 Th decay chain we report just the activity of the 
progenitor of the chain, because the chain is at secular 
equilibrium: in fact the crystal was grown more than 14 
years before the measurement, thus all the nuclides of the 
chain have the same activity. On the contrary for the 238 U 
chain we evaluate the activity of all isotopes since the 
equilibrium is broken. 



Chain 


Nuclide 


Activity 






[MBq/kg] 


232rp^ 


232 rp^ 


51±8 


238tj 


238tj 


< 10 




234 u 


< 7 




230 Th 


178±15 




226 Ra 


1403±43 




210p o 


(186±1)-10 3 


147 Sm 


147 Sm 


7±5 




1200 1400 1600 1800 2000 2200 2400 
Energy [keVJ 

Fig. 4. Low energy scatter plot. The green and blue lines 
represent the 3 a acceptance region, respectively, for /J/7 and 
a events. 



4.5 Half-lives of lead isotopes 

We report in Tab. [l] all the lead isotopes that we have 
investigated with our PbWC>4 crystal. The a//3 discrimi- 
nation allows us to easily disentagle the a and /J/7 particle 
interactions. Unfortunately, the more we go at low energy 
the less efficient is the discrimination power due to the 
poor energy resolution of the LD. For energy releases in 
the PbW04 crystal smaller than 1 MeV, /J/7 events leak 
in the a region, because the two bands start to merge. 
Thus /J/7 events induce a negligible background for the 
investigation of 204 Pb and 206 Pb decays, but not for 207 Pb 
and 208 Pb, which are at lower energy. This background is 
therefore evaluated by defining the 3 a acceptance region 
for a and /J/7 events, as shown in Fig. [4] 



4.5.1 



201 



Pb 



It is theorized that 204 Pb can a decay on 200 Hg, the Q- 
value of the transition is 1969.5 keV. The alpha energy 
spectrum in this energy region, shown in Fig. [5] clearly 
shows no peak, but we have a flat background, as al- 
ready mentioned in Sect. |4.4| We study an interval of 2.8 a 
(that corresponds to a confidence interval of 99.5%) cen- 
tered around the Q-value of the transition, being a the 
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1000 1200 1400 1600 1800 2000 



2200 2400 

Energy [keV ] 



Fig. 5. Energy spectrum of just a events, the dashed line 
represents the expected 210 Po flat background. The ROI for 
204 Pb and 206 Pb decays are highlighted. 



energy resolution of the 180 W peak (5.9 keV). The com- 
puted 210 Po background around the 2.8 a region of inter- 
est (ROI) of 16.4 keV, is 0.3 counts, while the observed 
number of counts is 0. Applying the Feldman-Cousins 
method [37] we are able to set a 90% C.L. limit on the 
half-life of 204 Pb, using the following formula: 



T\/2 — N nuc i e i 



N r , 



■ £tot , 



ln(2) (4) 



where: N nuc i ei is the number of nuclei of the analyzed 
isotope, N counts is the number of events estimated using 
the Feldman-Cousin approach, T meas is the measurement 
time in year and €tot is the total efficiency of event-based 
cuts, which is assumed to have the same value of the 180 W 
one. The computed lower- limit is: 



T? /2 ( 204 Pb) 



> 1.4 • 10 2U y at 90% C.L. 



0.5:- 



320 



340 



360 



380 



400 



420 



T? /2 r r Pb) 



> 1.9 



10 21 y 



at 90% C.L. 



4.5.4 208 Pb 



440 460 
Energy [keV 



Fig. 6. Low energy scatter plot analyzed for investigating 
207 Pb a decay. The blue and green lines represent the 90% 
acceptance region, respectively, for a and /3/j events. The red 
dashed line is the energy cut applied for selecting the region 
with the highest statistical significance, while the shadowed 
region is our ROI. 



we should define an energy cut that maximizes the signal- 
to-noise ratio at such low energy. In Fig. [6] is shown the 
scatter plot around the Q-value: 392.3 keV, and the ap- 
plied cut on the detected light that maximizes the statisti- 
cal significance in that region. The estimated cut efficiency 
is 53%. 

In the analyzed 2.8 a interval (16.4 keV) around the Q- 
value the number of observed events is and the estimated 
background is 1 counts, therefore we can set a lower-limit 
on the 207 Pb half-life: 



4.5.2 206 Pb 

206 Pb is energetically allowed to decay through the a chan- 
nel on 202 Hg. The energy released in the decay is 1135.5 keV. 
Looking at Fig. [5] we can estimate the 206 Pb half-life ap- 
plying the same procedure explained in the previous sec- 
tion. The events distribution in the energy spectrum is 
still flat in the 1 MeV region, so no excess is observed at 
1135.5 keV. The evaluated background is 0.2 counts, and 
the observed number of events in the 2.8 a energy interval 
(16.4 keV) around the Q-value is counts. The lower- limit 
on the 206 Pb half-life is: 



The last natural occuring lead isotope to be analyzed is 
208 Pb. It can a decay to 204 Hg with an energy transition 
of 516.9 keV. Also for this nuclide we apply the same pro- 
cedure previously explained. The studied a band is shown 
in Fig. [7J the cut applied on the detected light has an 
efficiency of 64%. 

The interval centered around the Q-value of the tran- 
sition is wide 2.8 a (16.4 keV). No events are observed in 
the ROI, and the estimated background is count. The 
lower- limit on 207 Pb half-life is: 



Tf /2 ( 208 P6) > 2.6-10 21 y at 90% C.L. 



(5) 



T? /2 ( 206 Pb) > 2.5 • 10 21 y at 90% C.L. 



5 Conclusion 



4.5.3 207 Pb 

As already explained, for low energy a decays, where events 
from the /? band can give a contribution to the a region, 



We successfully tested a 454.1 g PbW04 crystal as a bolome- 
ter. The crystal was grown using low background ancient 
Roman lead, and its performances and internal contami- 
nations were presented. 
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Fig. 7. Low energy scatter plot for investigating 208 Pb a decay. 
The blue and green lines represent the 90% acceptance region, 
respectively, for a and j3/^ events. The red dashed line is the 
energy cut applied for selecting the region with the highest 
statistical significance, while the shadowed region is our ROI. 

The separation achieved on the scintillation light is a 
powerful tool that allows us to reject events for the 
investigation of low energy rare a decays. The LY of a 
and (3/ j has been calculated. 

New more stringent upper limits on the half-lives of 
four natural lead isotopes has been estimated with a sen- 
sitivity of the level of 10 20 -10 21 y. 
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